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We propose a simple mean-field ansatz to study phase transitions from a topological phase to a 
trivial phase. We probe the efficiency of this approach by considering the string-net model in the 
presence of a string tension for any anyon theory. Such a perturbation is known to be responsible 
for a deconfinement-confinement phase transition which is well described by the present variational 
setup. We argue that mean-field results become exact in the limit of large total quantum dimension. 
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I. INTRODUCTION 


The blend of quantum computation and of topolog¬ 
ical phases of matter' have led to the idea of topo¬ 
logical quantum computation^""^. In this field, the es¬ 
sential ingredient is the construction of physical sys¬ 
tems sustaining exotic excitations known as non-Abelian 
anyons (see Ref.^ for a review). Being genuinely nonlo¬ 
cal, these anyons allow for efficient storage and manipu¬ 
lation of quantum information. Indeed, topologically or¬ 
dered systems^ are stable under local perturbations*’ and 
hence protected against undesirable effects such as de¬ 
coherence. However, strong enough perturbations, may 
drive the system to a nontopological phase. In recent 
years, many works have been devoted to the study of this 
robustness in microscopic models. Such an issue is diffi¬ 
cult to address since one has to deal with two-dimensional 
interacting quantum systems and the complex nature of 
the anyonic quasiparticles prevents one from using stan¬ 
dard methods. 

The goal of the present work is to propose a simple 
approach that may be considered as a mean-field theory 
for topological phases. To this end, we introduce a vari¬ 
ational ansatz which can describe topological as well as 
non topological phases. By construction, it also matches 
the exact ground state in some limiting cases. Thus, 
it aims at qualitatively describing phase diagrams while 
being quantitatively acceptable. Most models hosting 
topological quantum order are built as a sum of local 
commuting projectors (toric code^, string nets^,...). In 
lattice gauge theories, one often interprets these projec¬ 
tors as operators measuring effective fluxes and charges. 
The topologically ordered ground state (vacuum) is then 
defined as the flux-free and charge-free state. Elemen¬ 
tary excitations are obtained by locally violating this 
constraint. In two dimensions, excitations are pointlike 
anyons related by strings and their energy does not de¬ 
pend on their relative position so that topological phases 
are also called deconfined phases. A natural way to de¬ 
stroy topological order consists in adding a string tension 
that will drive the system to a confined phase. The pro¬ 


totypical Hamiltonian of such a system can be written 

V p I 

where Qy (Bp) are projectors measuring charges (fluxes) 
on vertices (plaquettes) of a two-dimensional graph and 
where L; is an operator acting on links which induces 
a string tension. In the following, we consider a two- 
dimensional plane with open boundary conditions so that 
the ground state is unique in the thermodynamical limit. 
Assuming non-negative couplings, the ground state of H 
is readily written in two limiting cases. On one hand, 
in the trivial phase Jv = Jp = 0, the ground state is a 
(polarized) product state denoted by |0), where all links 
are in the same state. On the other hand, for Jj = 0, 
the ground state is proportional to UyQvUp Bp\0). The 
main idea of our construction is to find a simple varia¬ 
tional state that bridges the gap between these two ex¬ 
treme cases. 

In this paper, we focus on the string-net model in the 
honeycomb lattice since it allows one to study a wide 
variety of topological phases^’®. Interested readers that 
are not familiar with this model can find a detailed study 
of this variational approach in the simpler case of the 
toric code model in Appendix B. 

The string-net Hamiltonian' is a special case of Eq. (1) 
where the operator Bp favors the zero-flux configuration 
in plaquette p. Here, we only consider states without 
charge excitation so that the Hilbert space is spanned by 
all link configurations satisfying the so-called branching 
rules (stemming from the fusion rules of the considered 
anyon theory). We thus drop the — Jv Qr term in the 
Hamiltonian. For simplicity, we also restrict our discus¬ 
sion to the string tension term introduced in Ref.** which 
involves Li operators enforcing a zero flux in link I of the 
lattice. Operators Bp and Li commute except if link I 
belongs to plaquette p. 

II. ANSATZ STATE AND ITS BASIC 
PROPERTIES 

To describe the phase transition separating the topo¬ 
logical phase from the trivial phase, we introduce the 
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FIG. 1. (Color online) Energy landscape as a function of a 
and X = Ji/Jp for — 2. The green line shows the position 
of the absolute minimum a{x). In this case, the transition is 
found to be continuous (second order). 


FIG. 2. (Color online) Energy landscape as a function of 
a and x — Ji/Jp for = 100. The green (dotted) line 
shows the position of the absolute (local) minimum a{x). For 
> 2 , the transition is found to be discontinuous (first or¬ 
der). 


following single-parameter variational state: 

|a) =Af 11(1 + aZp)\0), (2) 

p 

where 0 ^ a ^ 1, and Zp = 2Bp — 1 is such that 
Zp = 1. The normalization constant Af depends on the 
total quantum dimension D of the theory considered, 
on a, and on the system size (see Appendix A). Once 
again, the physical insight underlying this ansatz is that 
|a = 0) = |0) is the exact ground state for Jp = 0, while 
|a = 1) oc rip-^plO) the exact ground state for Ji = 0. 
Thus, one can expect that it captures the physics, at least 
qualitatively, for nonvanishing couplings. 

Interestingly, the structure of |a) implies that for any 
set Vn of n plaquettes, one has 


(n = n (-B^> 

pev,^ pG(P„ 


(l + o)^ 

Z?^(l — a)^ + 4a 


, ( 3 ) 


where (Oja = (a\0\a) (see Appendix A for details). 
This factorization property reveals the mean-field char¬ 
acter of I a). In addition, for Abelian theories, |a) can be 
rewritten as a simple product state in the dual plaquette 
(flux) basis. For illustration, let us consider the simplest 
Abelian theory, i.e., Z 2 (D^ =2). As shown in Ref.^'*, 
for this theory, the string-net model with a string tension 
can be mapped onto the transverse-field Ising model on 
the triangular lattice by setting XpXpi = 2Li — 1, where 
p and p' are plaquettes sharing link 1. In this dual rep¬ 
resentation, degrees of freedom are defined on plaquettes 
(instead of links) and operators Xp and Zp are the usual 
Pauli matrices. One can then compute the following ex¬ 
pectation values in the link basis (see Appendix A) 

and in the plaquette basis 

{Zp)e= COS0, {XpXp,)e=sinH={Xp)e{Xp,)e. (5) 

Here, we set |d) = (8)p[cos(0/2)|t)p + sin(d/2)|4,)p] where 
|t)p and IDp are the eigenstates of Zp with eigenvalues 


-1-1 and —1. Clearly, expressions (4) and (5) coincide 
provided a = tan(0/2). 

The Ztv case can be treated similarly by mapping the 
model onto the transverse-field iV-state Potts model'*’ 
(other models with Z 2 and Z 3 topological order have 
also been treated in the same vein""''’). Although no 
such mapping is known for non-Abelian theories (be¬ 
cause of the existence of multiple fusion channels), the 
state I a) can still be considered as a mean-field ansatz 
because of the factorization property (3). In other words, 
the present approach generalizes the canonical mean-field 
treatment (performed in the dual basis) implemented for 
Abelian anyons, to non-Abelian theories. 


III. RESULTS 

For any theory with total quantum dimension D, one 
can compute the variational energy per plaquette 

f ^ Jpfpia) + 3Jifi{a) 

e(Q() =-’ 

5(a) 

where 

g(a) = D^ [D'^{l-af +4a\^ , 
fp(a) = D^(l + a)2 [U2(l - af + 4a] , 

/i(a) = £i®(I - a)^ -f 8i:>'*a(l - af 

+24D‘^a^{l - a)2 -f 16a*’(2 - a). 

Details of the calculations are given in Appendix A. Set¬ 
ting X = Ji/Jp, the study of e(a) indicates that the sys- 
tern undergoes a phase transition at Xc = ■ Indeed, 

the minimum of e is obtained for a_ = I if x ^ Xc and for 
a+ < I if X ^ Xc (see Figs. 1-2 for illustration). At the 
transition, one has a+(xc) = This transition is 

second order for = 2 only, and first order for > 2. 

Interestingly, all these variational results only depend 
on D. This is reminiscent of the intrinsically local char¬ 
acter of the ansatz that does not take into account sub¬ 
tle effects due to nontrivial braiding statistics. Within 


( 6 ) 

( 7 ) 

( 8 ) 

(9) 
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Z 2 

Z 3 

Fibonacci 

Ising 

D'^ 

2 

3 

3.618 

4 

Xc (mean field) 

0.1667 

0.2222 

0.2412 

0.25 

Xc (series) 

0.2097'® 

0.2466'® 

0.261'' 

0.267'* 


TABLE I. Position of the transition point for several theo¬ 
ries computed with the mean-field ansatz ( 2 ) and with series 
expansions. 


this mean-field approach (single-plaquette approxima¬ 
tion), two theories with the same total quantum dimen¬ 
sion D are thus treated on an equal footing. Neverthe¬ 
less, from high-order series expansions, we know that, for 
instance, Z 4 and Ising theories = 4) have different 
ground-state energies^"*. 

Consequently, it is natural to wonder how these pre¬ 
dictions compare with exact results. First, it is worth 
noting that, in the topological phase {x < Xc), the en¬ 
ergy is minimized for a = 1 which is the exact result for 
a: = 0. For a = 1, the ground-state energy reads 


e(a_) 3x 

Jp 


( 10 ) 


which matches the exact small-a; perturbative expansion 
up to order 1 but does not give higher-order corrections. 
Secondly, in the opposite (large-a:) limit, the variational 
energy per plaquette can be expanded in powers of 1 /x 
and reads, at order 4, 

e(a+) _ 11 1 D"^-I I -ZD^ + 2 

Jy X x'^ 6D^ x^ 3611 ® 

1 2D^ - IID^ -p 19H2 - 10 

432118 ■ ^ 

This expansion matches the exact large-a: series expan¬ 
sion up to order 3 but not beyond. Once again, this is 
due to the local character of the ansatz that does not 
capture quantum fluctuations beyond a single plaquette. 

Another important remark concerns the behavior of 
the so-called Wilson loop operators denoted for a 
contour C„ enclosing n plaquettes and a string of type 
s (see Appendix A). In the deconfined (confined) phase, 
the expectation value of is expected to scale as the 
perimeter (area) of Remarkably, the present mean- 
field approach displays this behavior since 




4a 

zi^l^a)^Ta 


( 12 ) 


where and ds are the Frobenius-Schur indicator and 
the quantum dimension of the string s, respectively 
(see Appendix A). In the topological phase, one has 
(W^^)a_ = iCsds for any Cn, which can be interpreted 
as a trivial perimeter law with an infinite characteris¬ 
tic length. By contrast, in the polarized phase, one has 
(kF^ )a+ = where the characteristic area A 

is readily obtained from Eq. (12). 
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FIG. 3. (Color online) Comparison of variational results (up¬ 
per boundaries of shaded planes) with low- (high-) field series 
expansions shown in full (dashed) lines for theories discussed 
in Table I. Bare series at highest available orders^®^*® are 
displayed. 


Let us now compare the mean-field predictions with 
existing results. As explained above, for the Zjv theory 
{D^ = N)^ the model Is equivalent to the Wstate Potts 
model In a transverse field on the triangular lattice. This 
model is known to display a second-order transition (Ising 
universality class) for N = 2, and a first-order transition 
for ^ 3 (see Ref.^'* for a review). Thus, the present 
mean-field treatment gives the correct order of the tran¬ 
sition. In Table I, we give the position of the transition 
point Xc obtained from series expansions and from the 
present mean-field ansatz. Quantitatively, the difference 
between the results of both approaches decreases as 
increases. For the Potts model, the mean-field theory 
is even known to be exact for large In this 

limit, one obtains a first-order transition at Xc = 1/3. 
Since e(a) only depends on Z), this large-ZI mean-field 
result is expected to hold for all theories. 

However, for non-Abelian theories with finite H, the 
situation is more complex. In two recent studies^^ds^ ^g_ 
ing series expansion and exact diagonalizations, it has 
been claimed that the phase transition for Fibonacci and 
Ising theories is second order but the present mean-field 
approach predicts a first-order transition > 2). Al¬ 
though none of these methods are exact, we strongly be¬ 
lieve that a (weakly) first-order scenario is correct. Apart 
from the mean-field result, this conclusion relies on two 
observations which have been overlooked. 

The first one relies on strong similarities of the ground- 
state energy series expansions between Z 3 , Fibonacci, 
and Ising theories (see Fig. 3). In particular, a jump 
in the first derivative of the ground-state energy per pla¬ 
quette dejdx is observed at the transition. This jump 
was considered as an artifact due to a finite-order series 
in Refs.'^'d®. As can be seen in Fig. 3, the magnitude of 
this jump is found to increase with D, in agreement with 
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the mean-field result which yields 


dx 


X—Xc 


de 3(£|2 - 2)2 

x=x+ ~ Z)2(£)2 _ !)• 


(13) 


In addition, for > 2, the (relative) height of the en¬ 
ergy barrier at x = Xc between the two minima a± and 
the local maximum a* reads 


e{a±) — e{a*) 
e{a±) 




- 1 -h 4D2 - 4 
4 ( 1)4 + £)2 _ 1 ) ■ 


(14) 


In the limit {D^ — 2) <C 1, this relative energy vanishes 
as (£)2 — 2)4, indicating a weakly first-order transition. 
Such a behavior qualitatively explains why the transition 
for Fibonacci and Ising theories has been considered as 
second order in Refs, 8, 

The second argument that corroborates this scenario is 
based on the emergence of bound states in the low-energy 
spectrum inside the topological phase and will be dis¬ 
cussed elsewhere. Let us simply mention that such bound 
states are necessary although not sufficient to induce a 
first-order transition and they are present for > 2. 

The mean-field approximation can also be used to an¬ 
alyze the same model but in the ladder geometry for 
which several exact results are known’’ ’22-2 4 ^ 
dimensional case, the variational energy is straightfor¬ 
wardly obtained from (6) by merely replacing J\ by Ji/3. 
As for the two-dimensional case, the ansatz (2) predicts a 
second-order transition for = 2 and a first-order tran¬ 
sition for 1)2 > 2 although, for the ladder, the transition 
is known to be first order only if > 42 ’. Thus, the 
ansatz fails at describing the nature of the transition for 
Z)2 ^ 4, as already known for the Potts model2'’. Inter¬ 
estingly, the position of the mean-field transition point 
^.i^adder _ £>^-1 gggg I (self-dual point) in the large-D 
limit, which is the exact result for any Z)24. 


IV. CONCLUSION 


non-Abelian models, although the principles for doing 
so have been laid down^'’. The technique exhibited in 
the present paper can be considered as a first step, even 
though the tensor-network structure has been bypassed. 
Let us emphasize that single-parameter TNS have al¬ 
ready been proposed for Z 2 models2''’2’’, but their single- 
line structure leads to qualitatively wrong results (first- 
order transition). To solve this problem, Gu et al. in¬ 
troduced multiparameter double-line tensors2^. It seems 
that the double-line structure (or triple-line structure for 
non-Abelian theories) is crucial since it encodes informa¬ 
tion about plaquettes that is necessary for an area law 
in the confined phase. Our single-parameter ansatz sup¬ 
ports this conclusion. 

Let us also stress that tensors must be chosen carefully, 
in order to allow for topological states'^ ^^^^4 Por the toric 
code in a parallel magnetic field, we have shown (see 
Appendix B for a detailed calculation) that the topolog¬ 
ical entropy’'^’’^*’ vanishes for a < I, i.e., in the polarized 
phase, but is equal to — log 2 U (which is equal to —I 
since Z) = 2) for a = 1, i.e., in the topological phase. We 
conjecture that the same relations hold for the string-net 
model, for any theory. We leave the calculation of the 
topological entropy, or of other measures"*^’^® for future 
works. 

The use of TNS would furthermore allow one to 
study other transitions. For instance, for the Fibonacci 
theory’^, the ground state for Jp = 0 and Ji < 0 is the 
state |1) where all links carry a string 1, namely, a Fi¬ 
bonacci anyon. The transition from the string-net ground 
state to this state could thus be studied with a variational 
state I a) =AAnp(i -|-a^p)|l). However, analytical cal¬ 
culations are much harder in this case, so that numerical 
TNS methods would be extremely valuable. 

Finally, let us mention that it would be interesting to 
describe excitations in a variational setting and thus to 
study dynamical properties in the model, as was done 
in Ref.®”. We hope the present work will trigger such 
studies. 


To conclude, we would like to give some possible routes 
to go beyond the present approach. In Refs.2®’28j the 
ground state of the string-net model without string ten¬ 
sion has been written as a tensor-network state (TNS), 
involving a triple-line structure (that reduces to a double¬ 
line structure for Abelian theories). Following the steps 
detailed in these works, the state |a) could be written in 
the same way. The parameter a would only change the 
values taken by the tensors. Since |a) already captures 
semiquantitatively the physics of the transition induced 
by string tension, it seems reasonable to assume that per¬ 
forming a minimization over all parameters of the tensors 
should give more precise results. The first-order nature 
of the phase transition for > 2 should furthermore be 
favorable to the obtention of accurate results. 

TNS have already been successfully used to study 
phase transitions in Abelian models®’’2^^29^ However, 
they have not yet been applied to the more challenging 
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Appendix A: String-net model with string tension 
1. Definitions 

The Hamiltonian of the string-net model with a string 
tension is given by 

iL = -Jp^Rp-(Al) 

p I 

where the string-tension operator Li is the projector 
onto the trivial state |0); on the link 1. The opera- 
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tor Bp enforcing trivial flux in plaquette p is written as 
Bp = ^ dsBp. Here, we introduce dg = Xgdg, which 
is the product of the Frobenius-Schur indicator Kg and 
of the quantum dimension dg of the string s. The total 
quantum dimension is defined as D = °P" 

erator Bp inserts a string s in the links of plaquette p 
as defined in Appendix C of Ref 7. Since Bp acts as the 
identity on states satisfying branching rules (to which 
we restrict ourselves), we shall single it out and write 
Bp = 1. We introduce the operator Zp = 2Bp — 1 that 
satisfies Zp = 1 since Bp is a projector. With these no¬ 
tations and noting that do = 1, one obtains 

"*■ £)2 X ! (^ 2 ) 

s/O 

2. Normalization of |a) 

We consider the variational state 

|a) =AA n(> + aZp)\Q), (A3) 

p 

where 0 ^ a ^ 1 is a variational parameter. The fully 
polarized state |0) is defined as |0) = ®i|0);. The first 
task is to compute the normalization constant Af. To 
this end, let us note that 

2a 

{1+aZpY = {l+a^){l+r]Zp), with 77 = --(A4) 

1 -I- 

Thus, denoting iVp the number of plaquettes and since 
all Bp commute with one another, we find 

1 = {a\a) = AA2(1 -p a2)Afp^o| H (^ + ^Zp) |0). (A5) 

p 

For simplicity and since we are interested in the thermo¬ 
dynamical limit, let us assume open boundary conditions. 
Then, the only contribution to (oinp(i -|- r]Zp) | 0 ) comes 
from the term proportional to 1 that arises when expand¬ 
ing Upi^ + r]Zp). Indeed, the action of a on |0), 
for a boundary plaquette, introduces nontrivial strings in 
the boundary links that cannot be compensated by any 
other operator Bp^p. Using Eq. (A2), it is then easy to 
get the normalization condition 

1 = (a|a) = A/'^(l-I-a^)^pe^p, (A6) 


with 

n2 _ 2 

(A7) 


3. Computation of {Bp}a 

Let us pick a particular plaquette p and compute 
{Bp)a = {oi\Bp\a). Since all Bp commute with one an¬ 
other, we get 

{Bp)^=Af^{l + a^)^-{Q\Bp\{{l+riZp,)\Q). (A8) 

p' 

From the definition of Zp, it is easy to derive the iden¬ 
tity Bp{l + r]Zp) = (1 -I- r])Bp. The prefactor of 1 in 
this term is Proceeding along the same lines as for 
the normalization of |a), and using the expression of Af 
stemming from Eq. (A6), we then find 

{Bp).= ^-^-^. (A9) 


4. Computation of (HpTJpja 

Let Vn be a set of n plaquettes. The same argument 
as above shows that all plaquettes of Vn will have a con¬ 
tribution i, while other plaquettes have a contribu¬ 
tion 1. As a consequence 

(n< = (i^))” = w;. (AIO) 

P&Vn ^ 


5. Computation of {Li)c, 

Let us Anally turn to the computation of 
{Li)a = {ct\Li\a), which is a little bit more involved. We 
denote pi and p 2 the plaquettes sharing link 1. Then Li 
commutes with all Zp operators, except those acting at 
plaquettes pi and p 2 - As a consequence 


{Li)^=N^{\ + a^)^-^-^{Q\{l+aZp,){l+aZpPLi{l+aZp,){l+aZpP {l+r^Zp)\0) (All) 

pApi,P 2 


As for the previous two computations, the only contri- term proportional to 1 after expanding the operators, 
button to the matrix element (0| • • • |0) comes from the 
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Consequently, we can already take into account the con¬ 
tribution of rip^^pi as well as 

the expression of Af, stemming from Eq. (A5), to write 

(A12) 

where 

\^} = il+aZpJil+aZp^)\0). (A13) 
Denoting 1 -|- aZp = j31 + 'yCp, with 

o 1 ~ ^ J ry T r>s 

P = l-a , l=j^, and Cp = 

s^O 

(AM) 

one gets 

1^) = [/32i +c'pj +7'Cp,Cp,] |0). (A15) 

Since L/ enforces a trivial (s = 0) flux in link I, and 
since Cp operators introduce non-trivial fluxes, only the 
first and third terms of |i/') contribute to the matrix ele¬ 
ment appearing in Eq. (A12). When acting with Cp^Cp^ 
on |0), the only way to obtain a trivial flux in link I is 
to take the same s in Cp^ and in Cp^. Thus, one gets 
all possible states with a loop s surrounding plaquette pi 
and a loop s surrounding p 2 - Since one requires the link 
I to be in the trivial s = 0 state, the weight of these 
states is equal to ds as can be found by using Eq. (2.23) 
in Ref.‘*°. As a consequence, we obtain 

{'tp\Li\tjj) -1), (A16) 

so that 

(«»=(7^)'^[/3‘ + 7'‘(0^-1)1- (A17) 

6. Computation of e(o) = {H)^/Np 

Finally, we can compute the variational energy per pla¬ 
quette 

e(«) = = -MBp)o. - (A18) 

where the factor of 3 comes from the fact that on a 
honeycomb lattice, the number of links is three times 
the number of plaquettes. Replacing {Bp)a and (L/)a by 
their expressions, and simplifying everything, one gets 
the energy per plaquette given in the main text. Note 
that for the ladder geometry’^ the link operator Li only 
acts on rungs. As there are as many rungs as plaquettes, 

the variational energy per plaquette for the ladder reads 
gladder(^) = -Jp{Bp)^ - Jy{L{) 


7. Computation of {Wc^)a 


For a contour C„ enclosing n plaquettes, the Wilson 
loop operator Wq^ inserts a string s along C„. In princi¬ 
ple, one should consider two distinct operators, depend¬ 
ing whether the string lies above or below the lattice. 
However, for our ansatz state, these operators have iden¬ 
tical expectation values so that we denote both of them 
as . This operator is given by where 

fype-s simple-string operator defined in Ref.^, 
and is nothing but a multi-plaquette version of B^. As 

commutes with all Bp operators, 

+ a')'^'=(0| n (1 + vZp)W^jO). 

(A19) 

Since Wc^|0) is the state with a string s along C„, the 
only non-zero contribution comes from (Np — n) opera¬ 
tors 1 for plaquettes outside C„, and from n operators Bp 
inside C„, annihilating the string s (where s is the dual 
string of s). Each of the n fusions of s and s gives a 
factor Ksjds, and the resulting contractible s loop gives 
a factor ds- As a result 

(A20) 

Simplifying this expression finally yields 

= ><sds . (A21) 


Let us mention that, as for a single plaquette, one can 
build the projector b 

inside C„. This operator has the following expectation 
value: 




1 

IB 


l-b(D2 



(A22) 


From this expression, it follows that (Wcri)a=i = 1 as 
expected. Furthermore, when n = 1, one can check that 
(bkc„=i)a = {Bp)a given in Eq. (A9). This result allows 
one to rewrite the expectation value of Wilson operators 
as: 

= ^sds . (A23) 


Appendix B: Toric code in a magnetic field 


1. Definitions 

The Hamiltonian of the toric code in a magnetic field 
reads 


H = -jY,Av-jY,Bp-Y,h-(Ti, (HI) 

V p I 
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where h = {hx,hy^hz) is a uniform magnetic field and 
(T; = af) are Pauli operators at link I of a square 

lattice. Furthermore, Ay = Bp = 

where v and p, respectively, denote vertices and plaque- 
ttes of the lattice (see Ref.'^®, and references therein for 
a detailed discussion of this model). These operators all 
commute with one another and = Bp = 1. Note that, 
with these definitions. Ay and Bp are not defined as pro¬ 
jectors. In the following, we shall only consider a system 
in the thermodynamical limit with open boundary con¬ 
ditions. 


The structure of this state is simple enough to allow for 
straightforward calculations of all quantities appearing in 
the Hamiltonian. 


a. Normalization 


Since = 1, one has {1 + aAyY = {l + a'^){l+r]Ay), 
where 


■q = 


2a 

1 + a'^ 


(B4) 


2. Ansatz and limiting cases 


For a finite-size system with vertices and open bound¬ 
ary conditions, the normalization condition thus reads 


For a given direction of the magnetic field h, following 
the mean-field prescription previously detailed for string 
nets, we introduce the following variational state: 

no -h ttAy) no + l3Bp)\h), (B2) 

V p 

where \h) denotes the state fully polarized in the field 
direction. When h = 0, the exact ground state is given 
by a = /3 = 1, whereas for J = 0, it is obtained for 
a = /3 = 0. Although the normalization constant Af is 
hard to compute for arbitrary |/i), it is possible to find 
exact expressions for some particular field directions. 

In the following, we focus on two simple directions: 
the parallel-field case where the field points in the z (or 
equivalently x) direction'*and the transverse-field 
case where it points in the y direction. In the former 
case, a second-order phase transition in the Ising uni¬ 
versality class is known to occur for hz/J ~ 0.328**“"*®, 
whereas in the latter case, a first-order transition occurs 
for hy/J = I"**. As we will see, the ansatz state \a, P) 
qualitatively captures these two very different behaviors. 


3. Parallel field 

For a start, we consider a field h = (0, 0, h) pointing 
along the z axis. When J vanishes, the ground state 
is liy) = namely, the polarized state where all 

spins point in the z direction. In the opposite limit 
where /i = 0, the system is in the topological (toric code) 
phase. The ground state is then an eigenstate of all Ay 
and Bp operators, with eigenvalues 1, that can be written 
A/Tlt, liy)- For h ^ 0, the Hamiltonian still com¬ 

mutes with all Bp operators and the problem can then be 
mapped onto an Ising lattice gauge theory on the square 
lattice*'. The ground state is an eigenstate of all Bp’s 
with eigenvalues 1 which enforces /3 = 1 in Eq. (B2). 
Thus, we consider the following simple ansatz state: 


1 = (a|a) = Af^il + ( 55 ) 

V 

The only contribution to (til nt)(l + arises from 

the term proportional to 1 (i.e., that does not involve 
any Ay operator), since an Ay operator for a boundary 
vertex flip boundary spins. These spin flips cannot be 
compensated by the action of other Ay operators. As a 
consequence, the state |a) is normalized if the following 
condition holds: 

1 = {a\a) = Af"^ {\ + . (B 6 ) 

b. Computation of {Bp}c, and {Ay)a 

Since all Bp and all Ay operators commute, it is easy 
to see that 

(Bp), = {a\Bp\a) = 1. (B7) 

The computation of (A„), = {a\Ay\a) is also straight¬ 
forward since 

(A,), = Af\l + (^|A, [](1 + qAy,m 

v' 

= (tr|(? 7 l-f A„) (1+77A„/)|fi). (B 8 ) 

v' 

Here, we used the normalization condition (B 6 ) and the 
fact that Ay = 1. As for the calculation of the norm, the 
only nonzero contribution comes from the term propor¬ 
tional to 1 , so that one gets 

{Ay)a = q- (B9) 

This result is independent of and is thus valid in the 
thermodynamical limit. This will be the case for all quan¬ 
tities discussed below. 

c. Computation of (n„6v„ 


|a) = \a (3 =1) = A/'TT(1 -t-aA^)!-!)-) (B3) Fet V„ be a set of n vertices. The same argument as 

above shows that all vertices of V„ have a contribution q 




while other vertices have a contribution 1 , so that 


g. Analysis of the variational energy and phase diagram 


- / CK 


(BIO) 


•uGV„ 


This factorization property illustrates the mean-field 
character of the variational state |a). 


d. Computation of (ui) a 

We now turn to the calculation of {crf)a = {a\af\a) at 
link I: 

{'^i)a = A/'^('fi| ]^(1 + OiAy)ai JJ(1 + (Bll) 

V V 

We denote Vi and the two vertices that share link 1. 
Then afAy- = —Ay .af for j = 1,2, while erf commutes 
with all other Ay operators. Using the trivial identity 
(1 -I- aAy.){l — aAy.) = (1 — a^)!, we obtain 

(af)„=AA2(^| n {l+vAy)il-a^m. (B12) 

V^Vi ,V2 

We thus get 


(^f)a 


1-a^ 
1 + 


2 


(B13) 


e. Computation of {o'f)a and (erf )□, 


The variational energy e{6) exactly has the form one 
would obtain by (i) noting that the Hamiltonian is dual 
to the transverse-field Ising model on the square lat¬ 
tice, thanks to the duality transformation Ay = /if and 
erf = with vi and V 2 being the two adjoining ver¬ 

tices to link 1; (ii) performing a mean-field treatment from 
the dual spin-1/2 variables /r, namely, by writing the vari¬ 
ational state as a product state \0) = satisfying 

(/if)e = cos6», (/if)e = sin6», and = sin^ 0. Of 

course, a direct mean-field treatment based on the orig¬ 
inal variables a cannot describe this transition since a 
product state is topologically trivial. 

The energy e{9) can be studied easily. It has 
a single minimum at 0 = 0, i.e., a = 1, when 
X = h/ J ^ Xc = ^/^- When x > Xy, a. single minimum 
is found for cos0 = Xy/x, thus for a < 1. This shows 
that there is a second-order quantum phase transition 
at X = Xy, between the low-held (x < Xy) topological 
phase, and the high-held {x > Xy) polarized phase. The 
mean-held approach is thus able to capture the qualita¬ 
tive features of the phase transition, since it is known 
that the transverse-held Ising model has a second-order 
quantum phase transition at Xy ~ 0.328. However, the 
position of the critical point is about 24% oh since we 
hnd Xy = 1/4. 

It is also interesting to note that, in the topological 
phase, e{x < Xy) = —J which agrees with the order I 
perturbative expansion in the low-held limit h/ J <C 1. 
In the polarized phase, one has 


For the sake of completeness, let us mention that 

(aDa = (af)„ = 0, (B14) 

which follows from the fact that erf and tr/ both hip a 
single spin and that this single spin hip cannot be com¬ 
pensated by any product of Ay operators. 

/. Computation of the energy per link e 


J 

e(a; > Xc) =-/i-- - —, (BI7) 

which agrees with series expansion up to order 2 in the 
high-held limit J/h ^ 1. As explained for string nets, 
this is due to the fact that the mean-held ansatz only 
captures quantum huctuations at a single-vertex level 
which is not sufficient to obtain the exact contributions 
at higher orders. 


On the square lattice, in the thermodynamical limit, 
the number of plaquettes N-p equals the number of ver¬ 
tices A/,, and this number is half the number of links N\ 
of the lattice. As a consequence, the variational energy 
per link e(a) = {a\H\a)/N\ can be written as follows if 
one gathers all previous results 

e(«) = + ■ (^15) 

Denoting rj = = cos0 and (^(^2 = sin0, one hnally 

obtains 

e{9) =—'^(cosO + 1) — hs\v? 9. (B16) 


h. Topological entropy 

A reliable way to detect topological order in a given 
quantum state is to compute the topological entropy * 

In the following, we show that state |a) has a nonvanish¬ 
ing topological entropy only if a = I, which is in agree¬ 
ment with the nature of the phases expected on both 
sides of the critical point Xy. 

As shown in Ref.'*®, the topological entropy can be ex¬ 
tracted from the computation of the Renyi entanglement 
entropy. We consider a system with open boundary con¬ 
ditions, and we split it into two subsystems C and V, 
where C is simply connected, as shown in Fig. 4 for an 
example where C has a square shape. We denote c, d, 
and n the numbers of vertices fully included in C, fully 
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FIG. 4. (Color online) Partition of the lattice between two 
subsets C and T), where C is simply connected. 


included in T), and belonging to both C and I?, respec¬ 
tively. The aim is to compute the Renyi entanglement 
entropy between C and V when the system is in state |a), 
namely, £2 = ~ log 2 [Tr(pc)]j where pc = Trx)\a){a\- In 
the thermodynamical limit, and for a domain C that be¬ 
comes bigger and bigger, £2 = Pn — 7 -I- • • • where • • • 
contains all terms that vanish as n —>■ 00 . The first term 
is nonuniversal, contrary to the second one iStopo = — 7 , 
which is the topological entropy. It can be shown that 
•Stopo = — log 2 D, where D is the total quantum dimen¬ 
sion of the model under consideration. We shall prove 
below that 

*S'topo(a = !) = -! and 5'topo(0 ^ a < 1) = 0, (B18) 

so that only the state |a = 1 ) has topological proper¬ 
ties, with the expected quantum dimension D = 2, since 
there are four kinds of Abelian particles in the toric code 
mo del 

In order to perform this calculation, we begin by 
rewriting |a) as follows: 

|q) = J|(l -bQA^)|a)c O \a)v, (B19) 


defined similarly with domain V. Then, state |a) can be 
expanded as follows: 

\a) = mI |q!)c O \a)v A^|a)c O A®|a)x, 

V v£d 

+oP ^ 0 A® A® |a)x)-b ■ ■ • 

-bQ!"A^i • • • A^^ |q:)c ® A® • • • A®^ |q:)x>^ > (B23) 

where in the last term, all vi, - ■ ■ ,Vn are distinct and 
belong to d. 

For open boundary conditions, the states |a)x), 

{A^\a)-D,v & d), {A® A®Ja)x,,wi ^ U2 G 5}, 

A® • • • A^J\a)x> are all different and form an orthonor¬ 
mal set of states. It is thus easy to take the partial 
trace needed to compute the reduced density matrix 
Pc = Tr'p|Q!)(a|. One gets 

Pc = A4^ ( |q!)c c (q:| + A^ |a)c c(q:|A^ -b • • • 

V ved 

+a2”A^^ • • • A^^ \a)c c {a\A^^, ’ ■ ■ , (B24) 

where we did not write as many terms as before to keep 
things as readable as possible. 

To compute Tr(p^), one needs to find the spec¬ 
trum of Pc- For this, one has to see that all 
states appearing in the expression of pc, namely, 
W)c, {A^vW)c,v & d), {A^^A^Ja)c,'Ci ^'C 2 G 9}, •••, 
A^^ • • • A);^|a)c are normed but not orthogonal to each 
other, because one has the following identity 

<■■•< = 11^- (B25) 

vec 


where d denotes the boundary of C and V, namely, the 
vertices belonging to both C and V, and where 

\o)c =Afc + aA„)|tr)c with |fr)c = G);ec|t)G (B20) 

vGC 

\a)x> =A/x) (1 -b aA^)|'f|-)-p with jlD-p = ( 8 )iGX>lt)i(-B 21 ) 

vGTi 


Indeed, let us define two complementary states |'0i)c and 
IV' 2 )cj as states from the set written above that have the 
overlap 

c('0i|V'2)c = c(a|A^^ • • • A^^|a)c = c(a| H 

vec 

= (B26) 


These are normalized states, which impose the conditions 
(l-ba^)'^ = 1 and ^"^(14-0:^)^^ = 1 , that can be found 
as before, since C and V have (at least) one boundary. 
Knowing that JV^{1 -b a^Y+<^+'^ = we can rewrite 

\a)=M\{{l+aAY\a)c®\a)v, (B 22 ) 

v^d 

with A4^(l + — 1 - 

For each vertex z; G 9, we write Ay = AyA^, where 
the operator A^ = rizG^nC spins that are 

belong to both v and C, and where A® = Y\.i^vr\V af is 


the last expression being obtained as Eq. (BIO). Thus, 
two complementary states are not orthogonal. On the 
contrary, two states that are not complementary are eas¬ 
ily seen to be orthogonal. As a consequence, the density 
matrix pc has a block diagonal structure. Each block is 
a 2 X 2 matrix involving two complementary states of the 
form 

p^A =^2 (a'^Y'ipi)cc{'ipi \ +a^^”"'’^IV’2)cc(V’2|) , 

(B27) 

where j is the number of A^ operators appearing in state 
IV^i)c = Ay^ ■ ■ ■ Ay, |a), with 0 < j ^ n/2. For the sake 
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of simplicity, we shall consider that n is an even number, 


as is the case in Fig. 4. When j < n/2, there are 


such matrices. When j = n/2, there are 


matrices 



The matrices can be rewritten in an orthonormal 
basis made of states \(t>i)c and \ 4 ' 2 )c- For this, one first 
computes the overlap matrix O with matrix elements 
Ok,i = c(V’fc|'/’i)c with k and I taking values 1 or 2, 
namely. 


O = 




(B28) 


This symmetric matrix can be diagonalized by perform¬ 
ing a rotation, O = ^VW, with 


r = 


1 

71 



1 


and 


V = 


1 - 77 ° 

0 


0 

l + p" 


(B29) 

With these definitions, and noting that the diagonal ma¬ 
trix V has nonnegative diagonal elements, one can per¬ 
form the change of basis 


/|V’l)c\ ^t-ppl/2-p 7'/l)c\ _ 
\\i’ 2 )c j \\ 4 > 2 )c j 


One can then express p^^ in the orthonormal basis of 
states and \(j) 2 )c (the expressions being quite large, 
we shall not give them here). 

A check of the validity of the obtained expression is to 
compute Trpc. We find that 



domain C, c grows like n^). One then gets the following 
behavior of £2 ■ 


£2 


nlog2 


1 + 


(B34) 


where • • • represents terms that vanish when taking the 
limit n —>■ 00 . As a consequence, the topological entropy 
vanishes when 0 ^ a < 1. 


i. Wilson loops 

As already mentioned, the topological phase {x < Xc) 
is the deconfined phase of the Ising lattice gauge model, 
in which Wilson loops are known to obey a perimeter 
law^^b By contrast, in the polarized (deconfined) phase, 
these loops obey an area law. In the toric code model, 
Wilson loop operators Wc can be chosen as a product 
of erf operators along a closed contour, which is nothing 
but the product of all operators Ay surrounded by this 
contour. Using Eq. (BIO), one thus gets {Wc)a = 
where 77 = cos 9 = In the topological phase, the 

energy is minimized for a = 1 so that {Wc)a = 1 for any 
contour C. This can be interpreted as a trivial perimeter 
law with an infinite characteristic length. In the polarized 
phase {a < 1), one can write (Wc) = exp(—nln(l/77)), 
which is an area law with a characteristic area 1/ln(l/? 7 ). 
Our variational state thus correctly mimics the expected 
behavior of W in the deconfined phase as well as in the 
confined phase. 


4. Transverse field 

We now consider a field h = (0, h, 0) pointing along 
the y axis. In this case, the model is known to display 
a first-order quantum phase transition at the self-dual 
point h = ’ . The variational state reads 

|a,/3)=Afn(l + (y.Ay) IK* + /3iIp)K), (B35) 

V p 


as it should. Note that we have extended the sum over 
j = 0, • • • , n/2 to j = 0, • • • , n and have corrected the 
induced double counting by the prefactor 1/2. 

Similarly, one can compute Tr and deduce 


Tr {pD 





2a^ 
1 -I- 


n 


(B33) 



where |^) = (8)i|—>■);, is the polarized state where all spins 
point in the y direction. For the sake of completeness, 
we introduced two variational parameters a and (3 but, 
for symmetry reasons, we expect them to be equal. 

Since all calculations follow closely that of the preced¬ 
ing section, we shall directly give the results without fur¬ 
ther justification. The expectation values of charge and 
flux operators read 

{Av)a,p = ^ ^2 {Bp)a,l3 = ^2 ■ (B36) 


In the case a = 1, one finds Tr (p^) = 2^ ”, so that 
£ 2(0 = 1) = n — 1 and S'topo = —1, as expected'^®’^®. 
When 0 ^ a < 1, the second term in the above equation 
vanishes exponentially fast when n grows (for a generic 


In addition, one has {crf)a,p = 0, {(Tf)a,i 3 = 0, and 




f l-f 32 

\i + ay \l + P^ 


2 


(B37) 
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Setting cos6* = sin6* = and cost/) = 

T_o2 

sint/) = ijJ^ 2 , and keeping in mind that, in the thermo¬ 
dynamical limit, IVv = iVp = JV]/2, one gets the following 
energy per link: 

e(0, </)) = — -^ cos ^ ~ ^ cos 4> — h sin^ 9 sin^ (j). (B38) 

As expected, this variational energy is found to be min¬ 
imum for (j) = 0, i.e., for a = (3. Again, the expression 
of e{9,9) could have been obtained by using a duality 
transformation and treating the dual model in a mean- 
field way. The study of e{9, 9) shows that 0 = 0 is al¬ 
ways a minimum, but is the absolute minimum only for 
X = hj J < Xc = 27/32. For x ^ x* = a second local 


minimum appears and it becomes the absolute minimum 
for X > Xc- For x = Xc, two absolute minima coexist, at 
9 = 0 and at 0 = arccos(l/3). We thus find a first-order 
quantum phase transition aA, x = Xc- Consequently, our 
variational analysis is thus about 16% off the exact result 
Xc = I, and misses the self-duality of the model^'^. 

Finally, in the topological phase, e(x < Xc) = —J 
agrees with the low-field series expansion up to order 1 
in hj J. In the polarized phase, the series expansion of 
the variational energy at order 4 in J/h reads 

= (B39) 

which matches the high-held series expansion up to or¬ 
der 2 in J/h^^ (odd order contributions vanish). 
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